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Ch ' Abstract 

<tf ! 

. While the universe becomes more and more homogeneous at large 

scales, statistical analysis of galaxy catalogs have revealed a fractal 
structure at small-scales (A < 100/i -1 Mpc), with a fractal dimension 
7— i . D = 1.5 — 2 (Sylos Labini et al 1996). We study the thermodynamics 

of a self-gravitating system with the theory of critical phenomena and 

' finite-size scaling and show that gravity provides a dynamical mecha- 

■ nism to produce this fractal structure. We develop a field theoretical 
approach to compute the galaxy distribution, assuming them to be in 

qq | quasi- isothermal equilibrium. Only a limited, (although large), range 

ON ■ of scales is involved, between a short-distance cut-off below which 

^ cj other physics intervene, and a large-distance cut-off, where the ther- 

Q^-f modynamic equilibrium is not satisfied. The galaxy ensemble can be 

Q ■ considered at critical conditions, with large density fluctuations de- 

+3 , velopping at any scale. From the theory of critical phenomena, we 

^ | derive the two independent critical exponents v and r] and predict 

the fractal dimension D = 1/u to be either 1.585 or 2, depending on 
whether the long-range behaviour is governed by the Ising or the mean 
field fixed points, respectively. Both set of values are compatible with 
present observations. In addition, we predict the scaling behaviour 
of the gravitational potential to be r~2( 1+r ?). That is, r~ ' 5 for mean 
field or r~ 0,519 for the Ising fixed point. The theory allows to com- 
pute the three and higher density correlators without any assumption 
or Ansatz. We find that the iV-points density scales as r[ N 1 ^' D , 
when r\ >> ri, 2 < i < N. There are no free parameters in this 
theory. 

gravitation — galaxies : clusters: general — galaxies : general — galaxies 
: statistics — cosmology : miscellaneous — cosmology : large-scale structure 
of universe 
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1 Introduction 



One obvious feature of galaxy and cluster distributions in the sky is their 
hierarchical property: galaxies gather in groups, that are embedded in clus- 
ters, then in superclusters, and so on. (Shapley 1934, Abell 1958). Moreover, 
galaxies and clusters appear to obey scaling properties, such as the power-law 
of the two point-correlation function: 

£(r) oc r -7 

with the slope 7, the same for galaxies and clusters, of 1.7 (e.g. Peebles, 
1980, 1993). This scale-invariance has suggested very early the idea of fractal 
models for the clustering hierachy of galaxies (de Vaucouleurs 1960, 1970; 
Mandelbrot 1975). Since then, many authors have shown that a fractal 
distribution indeed reproduces quite well the aspect of galaxy catalogs, for 
example by simulating a fractal and observing it, as with a telescope (Scott, 
Shane & Swanson, 1954; Soneira & Peebles 1978). Sometimes the analysis 
has been done in terms of a multi-fractal medium (Balian & Schaeffer 1989, 
Castagnoli & Provenzale 1991, Martinez et al 1993, Dubrulle & Lachieze-Rey 
1994). 

There is some ambiguity in the definition of the two-point correlation 
function £(r) above, since it depends on the assumed scale beyond which 
the universe is homogeneous; indeed it includes a normalisation by the aver- 
age density of the universe, which, if the homogeneity scale is not reached, 
depends on the size of the galaxy sample. Once £(r) is defined, one can al- 
ways determine a length r where £(r ) =1 (Davis & Peebles 1983, Hamilton 
1993). For galaxies, the most frequently reported value is r ~ 5/iT 1 Mpc 
(where h = iJ /100km s _1 Mpc _1 ), but it has been shown to increase with 
the distance limits of galaxy catalogs (Davis et al 1988). r is called 'corre- 
lation length' in the galaxy literature. [The notion of correlation length £0 is 
usually different in physics, where £ characterizes the exponential decay of 
correlations (~ e~ r ^°). For power decaying correlations, it is said that the 
correlation length is infinite]. 

The same problem occurs for the two-point correlation function of galaxy 
clusters; the corresponding £(r) has the same power law as galaxies, their 
length r has been reported to be about r ~ 25h~ 1 Mpc, and their cor- 
relation amplitude is therefore about 15 times higher than that of galaxies 
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(Postman, Geller & Huchra 1986, Postman, Huchra & Geller 1992). The lat- 
ter is difficult to understand, unless there is a considerable difference between 
galaxies belonging to clusters and field galaxies (or morphological segrega- 
tion). The other obvious explanation is that the normalizing average density 
of the universe was then chosen lower. 

This statistical analysis of the galaxy catalogs has been criticized by 
Pietronero (1987), Einasto (1989) and Coleman & Pietronero (1992), who 
stress the unconfortable dependence of £(r) and of the length r upon the 
finite size of the catalogs, and on the a priori assumed value of the large- 
scale homogeneity cut-off. A way to circumvent these problems is to deal 
instead with the average density as a function of size (cf §2). It has been 
shown that the galaxy distribution behaves as a pure self-similar fractal over 
scales up to xs 100/iT 1 Mpc, the deepest scale to which the data are statisti- 
cally robust (Sylos Labini et al 1996; Sylos Labini & Pietronero 1996). This 
is more consistent with the observation of contrasted large-scale structures, 
such as superclusters, large voids or great walls of galaxies of ~ 200/?r 1 Mpc 
(de Lapparent et al 1986, Geller & Huchra 1989). After a proper statistical 
analysis of all available catalogs (CfA, SSRS, IRAS, APM, LEDA, etc.. for 
galaxies, and Abell and ACO for clusters) Pietronero et al (1997) state that 
the transition to homogeneity might not yet have been reached up to the 
deepest scales probed until now. At best, this point is quite controversial, 
and the large-scale homogeneity transition is not yet well known. 

Isotropy and homogeneity are expected at very large scales from the Cos- 
mological Principle (e.g. Peebles 1993). However, this does not imply lo- 
cal or mid-scale homogeneity (e.g. Mandelbrot 1982, Sylos Labini 1994): 
a fractal structure can be locally isotropic, but inhomogeneous. The main 
observational evidence in favor of the Cosmological Principle is the remark- 
able isotropy of the cosmic background radiation (e.g. Smoot et al 1992), 
that provides information about the Universe at the matter/radiation de- 
coupling. There must therefore exist a transition between the small-scale 
fractality to large-scale homogeneity. This transition is certainly smooth, 
and might correspond to the transition from linear perturbations to the non- 
linear gravitational collapse of structures. The present catalogs do not yet 
see the transition since they do not look up sufficiently back in time. It 
can be noticed that some recent surveys begin to see a different power-law 
behavior at large scales (A 200 — AOOh^ 1 Mpc, e.g. Lin et al 1996). 
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There are several approaches to understand non-linear clustering, and 
therefore the distribution of galaxies, in an infinite gravitating system. Nu- 
merical simulations have been widely used, in the hope to trace back from 
the observations the initial mass spectrum of fluctuations, and to test postu- 
lated cosmologies such as CDM and related variants (cf Ostriker 1993). This 
approach has not yet yielded definite results, especially since the physics of 
the multiple-phase universe is not well known. Also numerical limitations 
(restricted dynamical range due to the softening and limited volume) have 
often masked the expected self-similar behavior (Colombi et al 1996). A 
second approach, which should work essentially in the linear (or weakly non- 
linear) regime, is to solve the BBGKY hierarchy through closure assumptions 
(Davis & Peebles 1977; Balian & Schaeffer 1989). The main assumption is 
that the iV-points correlation functions are scale-invariant and behave as 
power-laws like is observed for the few-body correlation functions. Crucial 
to this approach is the determination of the void probability, which is a series 
expansion of the iV-points correlation functions (White 1979). The hierachi- 
cal solutions found in this frame agree well with the simulations, and with 
the fractal structure of the universe at small-scales (Balian & Schaeffer 1988). 
A third approach is the thermodynamics of gravitating systems, developped 
by Saslaw & Hamilton (1984), which assumes quasi thermodynamic equi- 
librium. The latter is justified at the small-scales of non-linear clustering, 
since the expansion time-scale is slow with respect to local relaxation times. 
Indeed the main effect of expansion is to subtract the mean gravitational 
field, which is negligible for structures of mean densities several orders of 
magnitude above average. The predictions of the thermodynamical theory 
have been successfully compared with N-body simulations (Itoh et al 1993), 
but a special physical parameter (the ratio of gravitational correlation energy 
to thermal energy) had to be adjusted for a better fit (Bouchet et al 1991, 
Sheth & Saslaw 1996, Saslaw k Fang 1996). 

We present in this article a new approach based on field theory and 
the renormalisation group to understand the clustering behaviour of a self- 
gravitating expanding universe. We also consider the thermodynamics prop- 
erties of the system, assuming quasi-equilibrium for the range of scales con- 
cerned with the non-linear regime and virialisation. We find an exact map- 
ping between the self-gravitating gas and a continuous field theory for a 
single scalar field with an exponential self-coupling. This allows us to use 
statistical field theory and the renormalisation group to determine the scal- 
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ing behaviour. The small-scale fractal universe can be considered critical 
with large density fluctuations developing at any scale. We derive the cor- 
responding critical exponents. They are very close to those measured on 
galaxy catalogs through statistical methods based on the average density 
as function of size; these methods reveal in particular a fractal dimension 
D ~ 1.5 — 2 (Di Nella et al 1996, Sylos Labini & Amendola 1996, Sylos 
Labini et al 1996). This fractal dimension is strikingly close to that ob- 
served for the interstellar medium or ISM (e.g. Larson 1981, Falgarone et al 
1991). We show in the present paper that the theoretical framework based 
on self-gravity that we have developped for the ISM (de Vega, Sanchez & 
Combes 1996a,b, hereafter dVSC) is also the dynamical mechanism leading 
to the small scale fractal structure of the universe. This theory is powerfully 
predictive without any free parameter. It allows to compute the iV-points 
density correlations without any extra assumption. 

We first clarify the definition of the average density we use for a fractal 
medium in §2, present the dynamical equations of the co moving fractal in 
§3, and apply our field-theory approach in §4. 



2 Correlation functions and mass density in 
a fractal 

The use of the two point correlation function £(r) widely spread in galaxy 
distributions studies, is based on the assumption that the Universe reaches 
homogeneity on a scale smaller than the sample size. It has been shown by 
Coleman, Pietronero & Sanders (1988) and Coleman & Pietronero (1992) 
that such an hypothesis could perturb significantly the results. The correla- 
tion function is defined as 

< n(ri).n(ri + r) > 

£( r ) = 5 1 

<n> 2 

where n(r) is the number density of galaxies, and < ... > is the volume 
average (over d 3 ri). The length r is defined by £(r ) = 1. The function £(r) 
has a power-law behaviour of slope —7 for r < r , then it turns down to zero 
rather quickly at the statitistical limit of the sample. This rapid fall leads 
to an over-estimate of the small-scale 7. Pietronero (1987) introduces the 
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conditional density 

_ < n(ri).n(ri + r) > 
~ <n> 

which is the average density around an occupied point. For a fractal medium, 
where the mass depends on the size as 

M(r) oc r D 

D being the fractal (Haussdorf) dimension, the conditional density behaves 
as 

r(r) oc r D - 3 

This is exactly the statistical analysis used for the interstellar clouds, since 
the ISM astronomers have not adopted from the start any large-scale homo- 
geneity assumption (cf Pfenniger & Combes 1994). 

The fact that for a fractal the correlation £(r) can be highly misleading 
is readily seen since 

T(r) 

C(r) = -M- ~ 1 
< n > 

and for a fractal structure the average density of the sample < n > is a 
decreasing function of the sample length scale. In the general use of £(r), 
< n > is taken for a constant, and we can see that 

D = 3 -7 . 

If for very small scales, both £(r) and T(r) have the same power-law be- 
haviour, with the same slope —7, then the slope appears to steepen for £(r) 
when approaching the length r . This explains why with a correct statistical 
analysis (Di Nella et al 1996, Sylos Labini & Amendola 1996, Sylos Labini 
et al 1996), the actual 7 ~ 1 — 1.5 is smaller than that obtained using £(r). 
This also explains why the amplitude of £(r) and ro increases with the sample 
size, and for clusters as well. 

In the following, we adopt the framework of the fractal medium that we 
used for the ISM (dVSC), and will not consider any longer £(r). 



3 Equations in the comoving frame 

Let us consider the universe in expansion with the characteristic scale factor 
a(t). For the sake of simplicity, we modelise the galaxies by points of equal 
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masses m, although they have a mass spectrum (it may be responsible for a 
multi-fractal structure, see Sylos Labini & Pietronero 1996). 

The present analysis can be generalised to galaxies of different masses 
following the lines of sec. IV in de Vega, Sanchez & Combes (1996b). We 
expect to come to this point in future work. 

If the physical coordinates of the particles are r, we can introduce the 
comoving coordinates x such that 

f— a(t) x 

The Lagrangian for a system of N particles interacting only by their self- 
gravity can be written as 



N 

i=l L 



—a(t) Xi 



m 
a(t) 



<l>(Xi(t)) 



(1) 



where 4>{x) is the gravitational potential in the comoving frame, determined 
by the Poisson equation 

V 2 (f)(x) = 4nG p(x, t) , (2) 
and p(x, t) is the mass density. For our system of point particles, 

N 

p(x,t) =m^25(x - Xi(t)) (3) 



and therefore the solution of the Poisson equation takes the form 



N 



<p(x) = —Gm 



i=l \ x x i{t)\ 

The canonical momenta and Hamiltonian of the system are 

Pi = m a{t) 2 Xi 
1 



H 



N 



n r 

E 

i=l L 



2ma(ty 



2ma(t) 2 

N 



Pi +777T<?W*)) 



N 2 Gm 2 



a(t) 

E 



fl W l<Kj<N 



(4) 



(5) 
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We see that the iV-particle Hamiltonian in cosmological spacetime eq.(^) 
can be obtained from the Minkowski Hamiltonian [a(t) = 1] by making the 
replacements 

m^ma(t) 2 , G^Ga(t)~ 5 . (6) 

As a first approximation, we shall assume in the following that the char- 
acteristic time of the particle motions under the gravitational self-interaction 
are shorter than the time variation of a(t). We can then consider that this 
system of self-gravitating particles is at any time in approximate thermal 
equilibrium. This hypothesis is true of course for structures that have already 
decoupled from the expansion, and are truly self-gravitating and virialised. 
It could be also valid for the whole non-linear regime of the gravitational 
collapse. As for the linear regime, we know already that the primordial fluc- 
tuations are not forgotten in the large-scale structures, and therefore the 
resulting correlations will depend on initial conditions, and not be entirely 
determined by self-gravity. 

The above assumption introduces a natural upper limit in the scales con- 
cerned by the theory developped below. Similarly to the case of the interstel- 
lar medium, the fractal structure considered is bounded by a short distance 
cut-off and by a large-scale limit as well (dVSC). 

The short distance cut-off corresponds to the appearence of other physics 
at short scale, essentially dissipative, which we do not need to introduce. 
In addition, the short distance cut-off avoids the gravo-thermal catastrophe. 
For the ISM, the cut-off was naturally the size of the smaller fragments, of 
the order of the Jeans length. Here the cut-off corresponds also to the size of 
the 'particles' considered, i.e. the galaxy size, below which another physics 
steps in, related to stellar formation and radiation. 

The present treatment can be generalized when thermal equilibrium only 
holds region by region (de Vega, Sanchez & Combes, in preparation). In such 
case we are lead to a quenched average over the temperature and we argue 
that the scaling properties are the same as in exact thermal equilibrium 
provided the temperature variations are smooth over the structure at the 
considered scale. 

The fact that in the catalogs, we are observing in projection large-scale 
structures at different epochs, with different values of the scale factor a(t), 
could slightly modify the fractal dimension. Even though fractal structures 
are self-similar, and scale-independent, the largest scales are systematically 
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observed at a younger epoch where the contrast has not grown up as high 
as today. This evolution effect however should be significant only at high 
redshift (> 1), and the present catalogs are not yet statistically robust so far 
back in time (the average redshift of optical catalogs is about 0.1). 

4 Application of renormalization group the- 
ory 

As in all scale-independent problems, where the fluctuations cannot be repre- 
sented by analytical functions, the renormalization group theory developped 
in the 1970's for the study of critical phenomena, appears here perfectly 
adapted (e.g. Wilson & Kogut 1974). We can consider the fractal structure 
of the Universe as the critical state of a system, where fluctuations develop 
at any scale, with a very large correlation length (asymptotically infinite). 
The fluctuations that are distributed as a fractal of dimension D are the 
large-scale structures of the universe (cf Totsuji & Kihara 1969). 

We have recently begun to tackle, with the tools of statistical field theory, 
the study of an iV-body system only interacting through their own self-gravity 
(dVSC). We have found an exact mathematical connection between this sta- 
tistical system and a scalar field with exponential self-coupling. We then 
used for it the powerful methods of field theory (e.g. Itzykson & Drouffe 
1989, Parisi 1988, Zinn- Justin 1989). Using the renormalisation group, the 
critical behaviour of this gravitational system has been described and its crit- 
ical exponents identified. This theory explains both the origin of the fractal 
structure, and predicts its fractal dimension D. This has been successfully 
applied to the interstellar medium (dVSC). Another approach has been pro- 
posed for galaxy correlations (not for ISM) (Hochberg & Perez Mercader 
1996), but it yields different critical exponents. 

Let us apply the theory to the system of galaxy points, already defined 
in the previous section. Since they are considered in approximate thermal 
equilibrium, we will use the grand canonical ensemble, that also allows a 
variable number of particles. The grand partition function of the system can 
be written as 
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oo Z N 



N=0 



N 

n 

1=1 



d 3 pi d 3 xi 
(2vr) 3 



-0H N 



(7) 



where z is the fugacity, which for an ideal gas of number density p is : 

/ , 2 \3/2 

z = Po ( 2-nmkT ) v 1 ls ^ ne Pl anc k constant). 

In dVSC we found a functional integral representation for the grand par- 
tition function 



Z 



J J V*e- 



■s[<K-)\ 



(8) 



i.e., Z can be written as the partition function for a single scalar field <fi(x) 
with local action 



S[cj>{.)\ 



1 



T. 



eff 



d x 



-(V0) 2 -/i 2 e<^ 



where 



-5/2 n m 2 

" zG(2m) 7 / 2 Vvr, T e// = 4vr 



In the 0-field representation, the mass density eq.(BI) is expressed as 



(9) 



(10) 



p(x) 



in 



T. 



eff 



T. 



eff 



(11) 



and the mass contained in a region of size R is 

/ 2 r R 



M{R) 



m ji 



eff 



& d 3 



X . 



(12) 



The mass parameter \x coincides at the tree level with the inverse of the Jeans 
length dj (dVSC) 



12 1 



(13) 



The functional representation for the grand partition function can be eas- 
ily generalized for an arbitrary scale factor a(t). After the changes specified 
above in eq.@, the local action becomes 



S[J>(.)} = 



a(t) 



T, 



eff 



d 3 x 



-(V0) 2 - ^a{tf 



(14) 
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Notice that all quantities depend on time through the scale factor a(t) only. 
There is no integration over t. 

The mass parameter \x in the 0-theory gets effectively multiplied by the 
scale factor a(t). Since the Jeans length dj ~ ^T 1 according to eq.(|13]), in 
comoving coordinates dj effectively becomes 

dj = \l — 



it fia(t) 

as one could have expected. 

On the other hand, the dimensionless coupling constant 

9 2 =^T eff 

is unchanged by the replacements of eq.@. 

Therefore, for any fixed time t we find the same scaling behaviour, after 
making the replacement 

li — > [A a(t) 
and keeping the coupling g unchanged. 



4.1 Scaling behaviour 

As is well known in the theory of critical phenomena (e.g. Wilson 1975, 
1983; Domb & Green 1976), physical quantities for infinite volume systems 
diverge at the critical point as A to a negative power, where A measures the 
distance to the critical point. The correlation length £o diverges as 

£o(A) ~ A~" , 

and the specific heat C behaves as 

C ~ A- a . (15) 

The critical exponents v and a are pure numbers that depend only on the 
universality class of the problem considered (e.g. Binney et al 1992). 

For a finite volume system, all physical quantities are finite at the critical 
point. And for a system whose size R is large, the physical quantities take 
large, but finite, values at the critical point. Thus, for large critical systems, 
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one can use asymptotically the infinite volume theory. In particular, the 

correlation length £ can be identified with the relevant physical scale R : 
£o ~ R- This implies that 

A ~ R~ l/v . (16) 

These concepts apply to the gravitational 0-theory since it exhibits scaling in 
a finite volume ~ R 3 (dVSC). Scaling behaviour was found for a continuum 
set of values of /i 2 and T e ff. 

We have previously shown (dVSC) that it is possible to identify 

A = 7^— = (271-mkT) 3 / 2 . (17) 

Notice that the critical point A = 0, corresponds to zero fugacity. Then, the 
partition function in the scaling regime can be written as 

Z(A)= f f P0 e -^+A/^^ , (i 8 ) 



where S* stands for the action at the critical point. 
We define the renormalized mass density as 

mp(x) ren = m e m (19) 

and we identify it with the energy density in the renormalization group (also 
called the 'thermal perturbation operator'). 

Since the 0-theory exhibits scaling (dVSC), the non-analytical part of the 
free energy is 

logZ(A) oc A 2 ~ a , 

so that its second derivative is C ~ A~ Q . Calculating the logarithmic deriva- 
tive of Z(h) with respect to A from eqs.(|T8|), using the standard relation 
between critical exponents in a three dimensional space a = 2 — 3v, and eqs. 



(|19"D and ([12]), we find that the mass fluctuations inside a volume of radius R 
(AM(R)) 2 = < M 2 > - < M > 2 , 

will scale as 

AM(R) ~ Ru . (20) 

The scaling exponent v can be identified with the inverse Haussdorf (frac- 
tal) dimension D of the system 

v 
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4.2 Critical exponents 



As usual in the theory of critical phenomena, there are only two independent 
critical exponents. All exponents can be expressed in terms of two of them: 
for instance the fractal dimension D = 1/u, and the independent exponent 
77, which usually governs the spin-spin correlation functions. The exponent 
r\ appears here in the 0-field correlator (dVSC), describing the gravitational 
potential, that scales as 



The values of the critical exponents depend on the fixed point that governs 
the long range behaviour of the system. 

The renormalization group approach applied to a single component scalar 
field in three space dimensions shows the presence of only two fixed points: 
the mean field point and the Ising fixed point. The scaling exponents as- 
sociated to the Ising fixed point are vi sing = 0.631..., D Ising = 1.585... , 
and rij sin g = 0.037... . The mean field value for the critical exponents are 

^meanf 2 ' ^meanf 2 and Tj mean j 0. 

The value of the dimensionless coupling constant g 2 = fjiT e ff should de- 
cide whether the fixed point chosen by the system is the mean field (weak 
coupling) or the Ising one (strong coupling). At the tree level, we estimate 
g ps -^L, where N is the number of points in a Jeans volume dPj. The coupling 
constant appears then of the order of 1, and we cannot settle this question 
without effective computations of the renormalisation group equations. At 
this point, the predicted fractal dimension D should be between 1.585 and 
2. 

4.3 Three point and higher correlations 

Our approach allows to compute higher order correlators without any extra 
assumption (de Vega, Sanchez & Combes, in preparation). 
The two and three point densities, 



< 4>(r) > ~ r 2 



± (!+*?) 



D(n,f 2 ) 



< n(fi) n(f 2 ) > 



£>(ri,f 2 ,f 3 ) 



< n{fi) n(f 2 ) n(f 3 ) > , 



(21) 
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can be expressed as follows in terms of the correlation functions: 

D{r 1 ,r 2 )=n l n 2 + C l2 (22) 
D(f 1 , r 2 , r 3 ) = n x n 2 n 3 + n x C 23 + n 3 C 12 + n 2 C 13 + C U3 . 

Here, 

rii =< n(fi) > , i = l,2,3, 

and Cij and are the two and three point correlation functions, respec- 
tively, 

Cij = C(fi, fj) 

The behaviour of n«, and C^k in the scaling regime follow from the 
renormalisation group equations at criticality (de Vega, Sanchez & Combes, 
in preparation). If we do not impose homogeneity at all scales, we find, 



< n(r) >~ A r 



D-3 



C(ri,r 2 ) ~ Brf , 

rt(-> -> - \ r i» r 2,»-3 „ 3(£)-3) /OQ \ 

G(ri,r 2 ,r 3 ) ~ Cr/ ; , (23) 

where A, B and C are constants and D = 1/v. 

We can now derive the three point density behaviour when one point, say 
fi, is far away from the other two. We find from eqs. fl22|) and (p3|), 



^ N r!»r 2 D _ 3 2(D-3) 

D(ri,r 2 ,r 3 ) ~ Ar± (n 2 n 3 + C 23 ) 

+B if D ~ 3) (n 2 + n 3 ) + C rf D -' i] (24) 

Notice that this expression is dominated by the first term since D — 3 < 0. 

Higher point distributions can be treated analogously in our approach. 
We find that the dominant behaviour in the iV-points density is 

C(f 1 ,f 2 ,...,r N )^ >>r ^ N r^ (25) 
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Notice that when homogeneity is assumed to hold over all scales, the critical 
behaviour of the iV-point correlation function involves a factor j-^ -1 ^- 0-3 ) ( 
(Itzykson & Drouffe, 1989). 

Eqs. (|24| - |25|) are qualitatively similar, although not identical, to the be- 
haviour inferred assuming the factorized hierarchical Ansatz (fhA), (Balian 
& Schaeffer 1989). That is, 



n 



~ n 



(l + br 
(l + rf 



D-3 
12 



D-3 



(26) 



D(ri,r 2 ,r 3 



,/fcA 



n 3 + n 3 b (r 



+n 6 Q 



n 



„D-3„D-3 
12 ' 13 



l + 6rg- 3 + 2rf- 3 (6 + g 3 Tf 3 



,D-3 
12 



D-3„D-3 
12 ' 23 



+ rf 3 3 4 / 



D-3 
23 



+ r 



D-3^,D-3 
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23 



D-3/ 



,D-3\ 



Qs^i 



2(D-3) 



where ria = |fi — f^l an d so on. b and Q% are constants. Notice that in the 
factorized hierarchical Ansatz, the fractal dimension D is not predicted but 
it is a free parameter. 

We see that the dominant behaviours in eqs.p^) and (|26|) are similar in 
case the scaling exponents D — 3 are the same. 



5 Conclusions 

The statistical analysis of the most recent galaxy catalogs, without the as- 
sumption of homogeneity at a scale smaller than the catalog depth, has 
determined that the universe has a fractal structure at least up to ~ 100/iT 1 
Mpc (Sylos Labini et al 1996). The analysis in terms of conditional density 
has revealed that the fractal dimension is between D = 1.5 and 2 (Di Nella 
et al 1996, Sylos Labini & Amendola 1996). We apply a theory that we 
have developped to explain the fractal structure of the interstellar medium 
(dVSC), which has the same dimension D. The physics is based on the self- 
gravitating interaction of an ensemble of particles, over scales limited both 
at short and large distances. The short- distance cut-off is brought by other 
physical processes including dissipation. The long-range limit is fixed by 
the expansion time-scale. In-between, the system is assumed in approximate 
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thermal equilibrium. The dynamical range of scales involved in this thermo- 
dynamic quasi-equilibrium is at present limited to 3-4 orders of magnitude, 
but will increase with time. 

The critical exponents found in the theory do not depend on the con- 
ditions at the cut-off, which determine only the amplitudes. The theory is 
based on the statistical study of the gravitational field: it is shown that 
the partition function of the N-body ensemble is equivalent to the partition 
function of a single scalar field, with a local action. This allows to use field 
theory methods and the renormalisation group to find the scaling behaviour. 
We find scaling behaviour for a full range of temperatures and couplings. 
The theory then predicts for the system a fractal dimension D = 1.585 for 
the Ising fixed point, or D = 2 in the case of the mean- field fixed point. 
Both are compatible with the available observations. The Appoints den- 
sity correlators are predicted to scale with exponent (N — 1)(D — 3) when 
ri »n, 2<i< N. That is, -(N-l) for the mean field, or -1.415 (iV-1) 
for the Ising point. 

We predict in addition a critical exponent — ^(1 + 77) for the gravitational 
potential: that is, —0.500 for mean field or —0.519 for the Ising fixed point. 
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